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Identifying key players in collective dynamics remains a challenge affecting a great variety of 
research fields, from the efficient dissemination of ideas to drug target discovery in biomedical 
problems. The difficulty lies at several levels: how to single out the role of individual elements in such 
intermingled systems, or which is the best way to quantify their importance. Centrality measures 
aim at capturing the influence of a node from its position in a network. The key issue obviated is 
that the contribution of a node to the collective behavior is not uniquely determined by the structure 
of the system but it is a result of the interplay between dynamics and network structure. Here we 
show that dynamical influence measures explicitly how strongly a node's dynamical state affects 
collective behavior. For diffusive processes in complex networks, dynamical influence quantifies 
how efficiently real systems may be driven by manipulating the state of single nodes. For critical 
spreading, it targets nodes with superior spreading capabilities. 



Complex networks are a groundbreaking concept that 
is helping to understand the behaviour of many chemi- 
cal, biological, social and technical systems [TJ [2] . Net- 
work representations are particularly suitable for systems 
where heterogeneity dominates and is crucial for dynam- 
ics [3J, where a few nodes are usually considered as the 
most important. Oftentimes, node importance is corre- 
lated with centrality measures, local [H [5] or global [S], 
which are generally based on a purely topological per- 
spective and do not explicitly account for the dynamics. 
However, dynamics is fundamental in assessing the im- 
pact of individual elements in global performance. Here, 
we show that dynamical influence is a centrality measure 
able to quantify how strongly a node's dynamical state 
affects the collective behavior of a system, taking explic- 
itly into account the interplay between structure and dy- 
namics in complex networks. We prove that it applies 
equally well to a variety of families of dynamical mod- 
els, from spreading phenomena at the critical point, for 
which dynamical influence targets nodes with superior 
spreading capabilities, to diffusive processes, for which it 
quantifies how efficiently real systems may be driven by 
manipulating the state of single nodes. 

Classical centrality measures in complex networks - 
like the degree or number of neighbors a node interacts 
with [5J [5] , betweenness centrality [7 counting the num- 
ber of shortest paths through a certain node, eigenvector 
centrality [5] based on the idea that relations with more 
influential neighbors confer greater importance, or the k- 
shell decomposition [9 that correlates with the outcome 
of supercritical spreading originating in specific nodes 
|10j- rely only on topology, even if an underlying process 
can be indirectly associated in some cases. In contrast, 
the impact of individual elements in the global perfor- 
mance of the system inevitably depends on the speci- 
ficities of the dynamics. Targeting individuals for vacci- 
nation strategies in epidemic processes is not the same 
as selecting electrical stimulation sites in the brain in 



order to suppress epileptic seizures. In this respect, a 
Laplacian-based centrality measure [PTUT3] . closely re- 
lated to PageRank [T3] . has been proposed recently to 
assess the importance of complex network nodes in spe- 
cific dynamical models. 

In this work, we provide a remarkably general and rig- 
orous framework where dynamical influence is defined 
as a centrality measure both on directed and on undi- 
rected complex networks and applies to a variety of fam- 
ilies of dynamical models, including epidemic spread- 
ing models like the susceptible- infected-removed (SIR), 
the susceptible-infected-susceptible (SIS), and the con- 
tact process, and diffusive processes like the voter model 
or phase coupled oscillators. In all cases, dynamical in- 
fluence is calculated as the left eigenvector associated to 
the largest eigenvalue zero of a characteristic matrix of 
the system, that encodes the interplay between topology 
and dynamics. 

More specifically, we focus on systems of N time- 
dependent real variables x — (x%, . . . ,Xjy) with coupled 
linear dynamics specified hy & N x N real matrix M 

x = Mx. (1) 

A first classification of the dynamics is obtained by con- 
sidering the largest eigenvalue ii max of M. For fi max < 0, 
x(t) converges to a null vector that represents a stable 
fixed point solution; for [i max > 0, indefinite growth from 
almost all initial conditions is observed. Suppose that M 
is such that a non-degenerate /z max = exists. Then, the 
scalar product <p c = c ■ x is a conserved quantity, where c 
is the left eigenvector of M for Umax, 

-Tf=c- = l cM \ ■ = ■ ( 2 ) 

The existence of the conserved quantity allows to cal- 
culate the final state in terms of the initial condition x(0) 
as 
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x(oo) :— lim x(t) 



x(0) 



c • e 



(3) 



where e is a right eigenvector of M for fi max . Eq. ([!]) 
implies that the projection of a;(0) on c is all the system 
remembers at large times about the initial condition x(0). 
Therefore, we call Ci the dynamical influence of element 
i in the dynamics under Eq. It is easily calculated 
without expensive numerical simulations. 

Let us first apply these insights to critical phenom- 
ena like spreading processes fTS] . In the SIR model [H] , 
each node is either susceptible, infected or removed. An 
infected node i transfers the epidemic along each of its 
outgoing arcs independently with probability /3; node i 
itself relaxes to the removed state at unit rate. As the 
spreading efficiency (details in Methods) of node i we 
define the expected fraction of nodes reached by an epi- 
demic outbreak initiated with a node i infected, all others 
susceptible. We call i the seed node. Viewing such a con- 
figuration as a small perturbation to the stationary state 
with all nodes susceptible, we approximate the dynamics 
by the linearization 



-x + f3A T x . 



(4) 



Here Xj (t) is the probability of node j to be infected at 
time t. The first term is the relaxation from the infected 
to the removed state at unit rate. The second term quan- 
tifies the transmission of the epidemic where the network 
enters by the transpose of its adjacency matrix A. 

Equation Q can be rewritten as Eq. ([I]) with M — 
f3A T — I and / the identity matrix. Matrix M has largest 
eigenvalue /i max — when the spreading probability j3 is 
the inverse of the largest eigenvalue of A, that is ft = 
fi c = l/a max . We take again c as a left eigenvector of 
M at \i max = or, equivalently, a right eigenvector of A 
at its maximum eigenvalue a max . Then the effect of an 
initial infection described by the probability vector x(0) 
is proportional to c • a;(0). 

Now we ask how well c may forecast the actual SIR 
spreading dynamics. Figure [T] shows that Ci is a good pre- 
dictor of spreading efficiency at critical parameter value 
ft = p c in a small social network. Dynamical influence Ci 
outperforms the predictions made by degree, shell index 
and betweenness centrality. Predictive power is quanti- 
fied by the rank order correlation (see Methods). 

Figure [2] shows the predictive power of dynamical in- 
fluence for spreading efficiency as a function of the in- 
fection probability in larger real-world networks and the 
Barabasi- Albert model. The results are as anticipated 
by the theory. Dynamical influence is a good predic- 
tor of spreading efficiency in the critical regime where 
(3/ fj c « 1. In this regime predictions by dynamical influ- 
ence outperform those by other quantities that are sup- 
posed to provide information about expected outbreak 
size. 

For sub- and supercritical parameter values [3, how- 
ever, the degree dj, of a node i is a good predictor of 
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FIG. 1: Spreading efficiency compared to centrality 
measures in a social network. The network of Zachary's 
karate club |17| has 77 edges connecting 34 nodes, here or- 
dered according to decreasing spreading efficiency. A mono- 
tonic decay of a centrality measure in the diagram indicates 
large predictive power for spreading efficiency. The rank order 
correlation of spreading efficiency is of 0.97 with dynamical 
influence, 0.86 with degree, 0.82 with shell index, and 0.79 
with betweenness centrality. Indexing of nodes is the same 
as in [17] . In the network drawing, circles and squares repre- 
sent the primary partitioning of the node set found by Girvan 
and Newman [18] . Spreading efficiency has been estimated at 
P = f3 c = 0.15 performing 10 6 independent runs of the SIR 
model per seed node. The largest eigenvalue of the network 
adjacency matrix is 6.65. 



spreading efficiency. In the subcritical regime, spreading 
is sparse and typically confined to the neighborhood of 
the seed node i, while in the supercritical regime, the 
epidemics rarely fails to spread to the whole system. In 
the critical regime in-between these extremes, infectious 
seeds are perturbations that trigger relaxation dynamics 
at all scales. This is reflected in a dynamics dominated 
by a marginal linear mode and a variety of possible fi- 
nal states. Dynamical predictions require then a global 
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FIG. 2: Predictive power for spreading efficiency as a func- 
tion of the spreading parameter (3 for the SIR model. Symbols 
are values of the rank order correlation coefficient of spread- 
ing efficiency with influence (squares), degree (triangles) and 
shell index (circles). Vertical dashed lines indicate the crit- 
ical value /3 C . The predictive power of betweenness central- 
ity (not shown) is below that of degree in all cases. The 
following networks have been used. E-mail interchanges be- 
tween employees of a university [19]; unweighted neural cir- 
cuitry of the roundworm C. elegans [201 121j ; snapshot of the 
Internet at Autonomous Systems level of Nov 08, 1997, see 
| http:/ /moat. nlanr.net | a realization of the Barabasi- Albert 
(BA) model of scale-free networks [4] with 1000 nodes and 
m = 2 edges added per node. Other realizations of the BA 
model yield qualitatively the same result. For the BA model, 
shell index is not a predictor because its value fej = m is the 
same for all nodes. For the neural network, being directed, 
out-degree instead of degree is used as a predictor and for 
calculating the shell index. 



view of the network structure and the final state is deter- 
mined by the conservation law associated with the lead- 
ing eigenvector c. The predictive power does not have 
not its maximum at the linear prediction f3/f3 c sa 1 pos- 
sibly due to the fact the infection cannot return to the 
node from which it has come |22) . 

Repeating the numerical experiments with the SIS 
model (see Methods) on the same networks, we obtain 
results (not shown) qualitatively similar to those of Fig- 
ure [2j Prediction by dynamical influence may also be 
applied to the contact process [23] • A is replaced by the 
stochastic adjacency matrix, the adjacency matrix after 
normalization such that each row sums up to 1. 

Coming back to the general framework Eq. [T] there 
is a class of dynamical processes in networks in which 
the property of M having a zero maximum eigenvalue 
appears naturally without the need of adjusting any pa- 




FIG. 3: Node dynamical influence and driving in a sys- 
tem of phase-coupled oscillators, a, Adaptation of the 
global phase of the four-node system in the inset when driving 
is applied to one of the oscillator nodes (an additional oscil- 
lator (not shown) coupled to the red circle node with fixed 
phase 7r/2. Adaptation is quick when driving at nodes 1 (o) 
and 4 (A), having high influence, and slow when driving at 
one of the other two nodes, having low influence, b, Driving 
efficiency (filled diamonds) and dynamical influence (shaded 
squares) for a system of phase oscillators coupled as the net- 
work of regions in the macaque cortex [24] . Driving efficiency 
of node i is measured as the time required to resynchronize 
with an additional input signal applied to a given node i. We 
say that the system has resynchronized when the global phase 
reaches ip(t) > (1 — e)x a with a tolerance e = 10 -2 , where the 
global phase %p(t) is computed as the argument of the global 
order parameter z(t) = r exp[iip(t)] = 2~Zi=i exp[iSi(t)]. For 
comparison, the degree ratio fcout/fcn of each node is also 
shown (open triangles). Since this is not an uncorrelated net- 
work, node influence deviates significantly from degree ratio. 



rameter. This is the case of diffusive processes denned 
by Eq. [TJwith M = —L and the Laplacian matrix entries 



-Ki 



N 

fc=l 



K„ 



ik 



(5) 



The zero eigenvalue of L is non-degenerate under mild as- 
sumptions |25| . For these processes our general analysis 
of Eq. [TJbecomcs exact. A prominent example of diffusive 
dynamics is the voter model |26| in which node i is in a 
spin state Sj e { — !,+!}. For this model, x; stands for 
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the ensemble average of spin i, Xi — (s,-), and gives 
the rate at which node % copies the state of node j. Differ- 
ent definitions of the voter model dynamics provide clear 
examples of how the concept of dynamical influence takes 
into account the interplay between topology and dynam- 
ics: For link update dynamics in an undirected network, 
an ordered pair of nodes is chosen in each step and 
node i copies the state of node j. The rate matrix K 
becomes then the adjacency matrix A. As a consequence 
Ci = 1/N, the average magnetization $D i=1 CiXi is con- 
served, and all nodes have the same dynamical influence 
independently of the topological features of the network. 
In the more standard node update voter dynamics, at each 
step one node i (having degree di) is selected at random 
and copies the state of one of its neighbors j, also selected 
at random. In this case Kij = Aij/di, so that Kij is no 
longer a symmetric matrix, the conserved quantity is a 
weighted magnetization 27] and the dynamical influence 
of node i is proportional to its degree di. 

For diffusive processes, the system is driven towards 
a homogeneous final state with x* :— Xi{oo) = Xj(po) 
for all i and j. Although x* takes continuous values, 
each realization of the voter dynamics in a finite system 
eventually reaches a homogeneous absorbing state with 
either all nodes in the state +1 or all in the state —1. The 
influence Ci of a node weights the initial state of node i 
in the probability of which of the two absorbing states 
is reached. When all nodes arc equivalent x* is just the 
average of the initial values of the nodes, but otherwise 
x* is given by the weighted average of the two absorbing 
states with the probability of reaching each of them. The 
value Ci has an alternative interpretation as a stationary 
density of a random walk |12j . 

The meaning of dynamical influence also manifests it- 
self in the practical task of driving a system efficiently. 
In the context of the voter model, this task might be 
phrased in terms of the zealot problem [25J . One consid- 
ers a special directed network in which a given node (the 
zealot) does not copy the state of any of its neighbors. 
The question is the efficiency of the zealot in driving all 
other nodes to the zealot state. To show the broad appli- 
cability of the dynamical influence concept, we address 
this question of driving efficiency considering the problem 
of phase-coupled oscillators described by the Kuramoto 
model [29 . Assuming all oscillators have the same in- 
trinsic frequency u, the phase variable Xi of oscillator i 
advances as 

JV 

Xi = u! + Kij sin(a;j — Xi) . (6) 
i=i 

with a matrix K of non-negative coupling strengths. 
Around the synchronized state, phase differences are 
small. Without losing generality, we choose lj = 0. By 
approximating each sin-term with its argument, a linear 
homogeneous system as in Eq. is recovered. 

We study a scenario with initially all oscillators i in 
phase Xi(0) = 0. An additional node a with constant 



phase x a — 7r/2 is added to the system and linked 
through an additional edge to a chosen node i. We mea- 
sure the time Tj the system takes to reach the new ho- 
mogeneous state with Si = ir/2 for all nodes i. The 
dynamical evolution of these systems is illustrated by 
studying the motif in the inset of Fig. [3^,. The global 
phase ip(t) converges faster to the external forcing when 
the driving is applied to the nodes with higher influence, 
and the convergence of the different nodes depends on 
their relative network position in relation to the driver. 
In Fig. |3]d, we show the results on a directed network 
of phase oscillators connected as the network of regions 
in the macaque cortex |24) . The extremely high correla- 
tion between dynamical influence and driving efficiency 
clearly shows that influence is an excellent proxy to iden- 
tify better targets for controlling global behavior, even 
in non-linear dynamical systems. The extremely high 
predictive power of dynamical influence for driving effi- 
ciency (rank order correlation of driving efficiency with 
dynamical influence is 0.97, while 0.66 with degree ratio 
and —0.09 with betweenness) clearly shows that influ- 
ence is an excellent proxy to identify better targets for 
controlling global behavior, even in non-linear dynamical 
systems. 

To summarize, dynamical influence is a general cen- 
trality measure which is not solely based on topology 
but specifically designed for dynamical processes on com- 
plex networks. For critical epidemic spreading, dynami- 
cal influence is a good predictor of spreading capabilities. 
For diffusion, dynamical influence quantifies the impact 
of the dynamical states of single nodes on the asymp- 
totic homogeneous state, and beyond that, it proves to 
be a high-quality proxy for driving efficiency, uncover- 
ing which are the best target nodes in real networks to 
be forced in order to drive the system towards specific 
states. In a broader context, the identification of these 
targets has fundamental implications and practical ap- 
plications on strategies with an interest in controlling 
collective behavior, from social influence to biomedical 
responses. 



METHODS 

Centrality measures 

The degree di of node i is the number of nodes i is 
connected to. In directed networks, in- and out-degree 
d\ n and d° ut are distinguished. For the matrix averag- 
ing over all adjacency matrices of networks with fixed 
node degrees, Ci = di is a left eigenvector for the largest 
eigenvalue. Likewise, the degree ratios d° ut /d\ a form a 
left eigenvector of the Laplacian matrix averaging over 
all networks with given node degrees [30l EI] . 

The betweenness centrality bi of a node i quantifies the 
fraction of shortest paths that pass through this node [7] . 
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It is defined as 



U;k) 



(7) 



where the summation runs over all ordered node pairs 
(j, fc); ajk denotes the total number of shortest paths 
from node j to node fc; (Tjk{i) is the number of such 
paths running through node i. 

The shell index ki [1011311 123] of a node i is derived from 
the consideration of the fc-core [33] for integer k > 0. The 
fc-core of a network is the largest induced subnetwork in 
which all nodes have degree at least k. Starting from 
the full network, the fc-core is obtained by deleting nodes 
(together with their edges) with degree strictly less than 
k until no such nodes are left. The shell-index ki is the 
largest value k such that node i is contained in the fc- 
core. In case of directed network, the fc-core is based on 
the out-degree. 



Epidemic models 

We simulate the SIR model of epidemic spreading in 
the time-discrete version. Transitions between the three 
states (S,I,R) are as follows. If node i is in the S (sus- 
ceptible) state and has v infected (I) neighbors at time t, 
then node i remains susceptible with probability (1 — (3)" , 
otherwise i is infected at time t + 1. If node i is in the 
infected state at time i then i is in the R (removed) state 
at time t + 1. In the SIS model, at difference with SIR, a 
node infected at time t is susceptible again at time t + 1. 
The probability of being removed in the SIR model does 
not enter in the linearized Equation (4| because it ap- 
pears only in a second order term in the equation for x. 
Therefore Equation Q gives the same linear description 
for the SIR and SIS models. 



The system is in an absorbing configuration if none of 
the nodes is infected. For both models, outbreak size is 
the number of nodes having been infected at least once 
before reaching an absorbing configuration. The spread- 
ing efficiency of node i is the average outbreak size when 
initiating the dynamics with node i infected and all oth- 
ers susceptible. 



Rank order correlation 



by 



For a vector x € M. n , the rank of component i is given 



n(x) = 1 + \{j ^ i\xj > Xi}\ + - \{j ^ i\xj 



*i}\ (8) 



The rank order correlation coefficient p(x, y) between two 
such vectors x and y is the Pearson correlation coefficient 
between the rank vectors r(x) and r(y). Thus p(x,y) 
takes values in [—1, 1] with p(x, y) = +1 ( — 1) if and only 
if x and y are in a strictly increasing (decreasing) relation. 



ACKNOWLEDGMENTS 

K.K. acknowledges financial support from Volkswa- 
genStiftung and from European Commission NEST 
Pathfinder initiative on Complexity through project 
SYNLET (Contract 043312). M.A.S. acknowledges fi- 
nancial support by the Ramon y Cajal program of 
the Spanish Ministry of Science, MICINN Project No. 
BFU2010-21847-C02-02, and Generalitat de Catalunya 
grant No. 2009SGR1055. V.M.E. and M.S.M. ac- 
knowledge financial support from MEC (Spain) through 
project FISICOS (FIS2007-60327). 



[1] S. N. Dorogovtsev and J. F. F. Mendes, Evolution of 
networks: From biological nets to the Internet and WWW 
(Oxford University Press, Oxford, 2003). 

[2] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

[3] R. Pastor-Satorras and A. Vespignani, Physical Review 
Letters 86, 32003203 (2001). 

[4] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[5] H. Jeong, S. P. Mason, A.-L. Barabasi, and Z. N. Oltvai, 
Nature 411, 41 (2001). 

[6] S. Wuchty and P. F. Stadler, Journal of Theoretical Bi- 
ology 223, 45 (2003). 

[7] L. Freeman, Sociometry 40, 35 (1977). 

[8] P. Bonacich, Journal of Mathematical Sociology 2, 113 
(1972). 

[9] B. Bollobas, Graph Theory and Combinatorics: Pro- 
ceedings of the Cambridge Combinatorial Conference in 
Honor of P. Erds Vol. 35 (Academic, 1984). 
[10] M. Kitsak, L. K. Gallos, F. Liljeros, L. Muchnik, H. E. 



Stanley, and H. A. Makse, Nature Physics 6, 888893 
(2010). 

[11] N. Masuda, Y. Kawamura, and H. Kori, New Journal of 
Physics 11, 113002 (2009). 

[12] N. Masuda and H. Kori, Phys. Rev. E 82, 056107 (2010). 

[13] B. Golub and M. O. Jackson, American Economic Jour- 
nal: Microeconomics 2, 112 (2010). 

[14] S. Brin and L. Page, Proceedings of the Seventh Interna- 
tional World Wide Web Conference, Computer Networks 
and ISDN Systems 30, 107 (1998). 

[15] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, 
Rev. Mod. Phys. 80, 1275 (2008). 

[16] J. D. Murray (2004). 

[17] W. W. Zachary, Journal of Anthropological Research 33, 
452 (1977). 

[18] M. Girvan and M. E. J. Newman, Proceedings of the Na- 
tional Academy of Sciences of the USA 99, 7821 (2002). 

[19] R. Guimera, L. Danon, A. Diaz-Guilera, F. Giralt, and 
A. Arenas, Phys. Rev. E 68, 065103 (2003). 



6 



[20] D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). 

[21] J. G. White, E. Southgate, J. N. Thomson, and S. Bren- 
ner, Philosophical Transactions of the Royal Society of 
London. B, Biological Sciences 314, 1 (1986). 

[22] R. Parshani, S. Carmi, and S. Havlin, Phys. Rev. Lett. 
104, 258701 (2010). 

[23] T. E. Harris, Ann. Probab. 2, 969 (1974). 

[24] C. J. Honey, R. Kotter, M. Breakspear, and O. Sporns, 
Proceedings of the National Academy of Sciences of the 
USA 104, 10240 (2007). 

[25] R. Agaev and P. Chebotarev, Linear Algebra and its Ap- 
plications 399, 157 (2005), special Issue devoted to pa- 
pers presented at the International Meeting on Matrix 
Analysis and Applications, Ft. Lauderdale, FL, 14-16 De- 
cember 2003. 

[26] R. Holley and T. Liggett, Annals of Probability 3, 643 



(1975). 

[27] K. Suchecki, V. M. Egufluz, and M. San Miguel, Euro- 
physics Letters 69, 228 (2005). 

[28] M. Mobilia, Phys. Rev. Lett. 91, 028701 (2003). 

[29] J. A. Acebron, L. L. Bonilla, C. J. Perez Vicente, F. Ri- 
tort, and R. Spigler, Reviews of Modern Physics 77, 137 
(2005). 

[30] M. A. Serrano, K. Klemm, F. Vazquez, V. M. Egufluz, 
and M. San Miguel, Journal of Statistical Mechanics: 
Theory and Experiment 2009, P10024 (2009). 

[31] N. Masuda and H. Ohtsuki, New Journal of Physics 11, 
033012 (2009). 

[32] S. B. Seidman, Social Networks 5, 269 (1983). 

[33] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, 
Phys. Rev. Lett. 96, 040601 (2006). 



